Quantum correlated light pulses from sequential superradiance of a condensate 
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We discover an inherent mechanism for entanglement swap associated with sequential superradi- 
ance from an atomic Bose-Einstein condensate. Based on careful examinations with both analytical 
and numerical approaches, we conclude that as a result of the swap mechanism, Einstein-Podolsky- 
Rosen (EPR)-type quantum correlations can be detected among the scattered light pulses. 
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I. INTRODUCTION 

Superradiance (SR) commonly refers to cooperative 
emission from an ensemble of excited atoms with initial 
coherence or from an ensemble of radiators with an ini- 
tial macroscopic dipole moment. As coherence enhanced 
radiation, SR was introduced by Dicke Jjli] in 1954, and 
first observed experimentally in 1973 It occurs in 
many systems from thermal gases of excited atoms 
4 1 and molecules 0, quantum dots and quantum wires 
E Hi 1 to atomic Bose-Einstein condensates (BEC) [1] , 
Rydberg gases Q, and molecular nanomagnets [13] ■ Re- 
cently, serious efforts have been directed toward the study 
of quantum entanglement between condensed atoms and 
SR Ught pulses [llj and entanglement between atoms 
through SR Several promising appHcations, including 
prospect for quantum teleportation in entangled quan- 
tum dots via SR, are proposed 6']. 

In the pioneering experiment of SR from an elongated 
condensate, a continuous wave (cw) pump laser intersects 
along the short transverse direction [8]. The scattered 
radiation is dominated by axial or the so-called end-fire 
modes fT*]. The atoms experience recoils as a result of 
the momentum conservation, exhibiting a fan-like pat- 
tern, which refiects the condensate side-mode distribu- 
tion. More recently, the Kapitza-Dirac regime of SR was 
observed [T^ in a pulsed pump scheme, with momentum 
side-modes displaying the characteristic X-shaped pat- 
terns. In this regime, it is predicted that SR pulses must 
contain quantum entangled counter-propagating photons 
from the end-fire modes [isj . It was proposed that quan- 
tum entanglement arises from correlations of backward 
and forward scattered atoms and from the interplay be- 
tween optical and atomic fields [iS^. In this study we 
show that even for a cw pumped condensate with scat- 
tered atoms forming a forward fan-like pattern, quantum 
entanglement of the end-fire modes still exists due to an 
entanglement swap mechanism which we clearly identify 
during sequential SR process. In quantum information 
language, entanglement swap is a techni que to entangle 
particles that never before interacted p^. llSl. [TgI . \vf\ . 

Sequential SR involves successive scattering of the 
pump laser from the initial momentum distribution of 



a condensate [8| . Previous studies on SR from an atomic 
gas have observed multiple pulses or ringing effects, es- 
pecially among dense atomic samples. Ringing is often 
explained in terms of the pulse propagation effect p^ . 
where the finite size and shape of the medium plays sig- 
nificant roles [l^, HO]. Adopting semi-classical theories, 
detailed modeling of SR from atomic condensates have 
been very successful, essentially capable of explaining 
both spatial and temporal evolutions of atomic and opti- 
cal fields [21I, m, [H, 111] . The semi-classical treatments, 
however, can account neither for the influence on se- 
quential scattering associated with ring from side-mode 
patterns nor for quantum correlations between end-fire 
modes. 

In this paper, we investigate Einstein-Podolsky-Rosen 
(EPR)-type [2^, [2^ quantum correlations between end- 
fire modes. Such correlations can be detected with well 
known methods developed for continuous variable entan- 
glement in down-converted two- photon systems [13, [l^ , 
employing equivalent momentum and position quadra- 
ture variables as observable. 

The paper is organized as follows. In sec. |TT1 we in- 
troduce the relevant concepts and describe the model 
system we consider for investigating sequential SR. We 
identify the various approximations and derive the full 
second quantized effective Hamiltonian. In sec. IIIIl we 
review the criteria for continuous variable entanglement, 
with which we confirm the existence of quantum corre- 
lation between SR photons from the end-fire modes. In 
sec. IIV[ we analytically solve the effective Hamiltonian 
under parametric and steady state approximations. We 
clearly identify the swap mechanism, and intuitively ex- 
plain the steps involved for the model Hamiltonian to 
generate EPR pairs out of non-interacting photons. This 
represents the key result for this article. In sec. |Vl we 
describe the method of our numerical calculations under 
a proper decorrelation approximation. The results are 
presented in sec. IVIl where we first examine the tempo- 
ral dynamics of the entanglement in connection with the 
accompanying field and atomic populations. This helps 
to illustrate the swap of atom-photon entanglement to 
the photon-photon entanglement. We then study care- 
fully this swap effect, introduce the effect of decoher- 
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ence, and consider the effect of SR initialization from a 
two-mode squeezed vacuum and the dependence on the 
increase/decrease of number of atoms. Sec. IVIII contains 
our conclusion. 



II. SEQUENTIAL SUPERRADIANCE AND THE 
EFFECTIVE HAMILTONIAN 

In this section, we briefly describe the unique prop- 
erties of SR, e.5., the directional and sequential nature 
of the emitted pulses. We will introduce the concept of 
sequential SR, in terms of what occurs in an elongated 
cigar-shaped atomic condensate. We derive the second 
quantized effective Hamiltonian, where the optical fields 
are treated quantum mechanically, in order to take into 
account the interaction of all side-modes with a common 
photonic fields. 




Laser Beam 



FIG. 1: (Color online) A fan-like atomic side-mode pattern 
up to second order sequential superradiant scattering. 



A. Sequential SR 

We consider an elongated condensate, of length L and 
width W , that is axially symmetric with respect to the 
long direction of the z-axis. It is optically excited with 
a strong pump laser of frequency cjq, detuned from the 
atomic resonance frequency loa by A = oja — wo. The 
laser beam is directed along the y-axis, perpendicular to 
the long axis of the condensate, and linearly polarized in 
the a;-axis. 

When the pump laser is sufficiently strong, the oc- 
cupation of atoms in the excited state becomes macro- 
scopic, beyond the threshold for collective emission. The 
excited atoms, interacting through the common electro- 
magnetic field, start to make collective spontaneous emis- 
sion pj]. During this dynamic process, the uncorrelated 
state of our model system first evolves linearly into a 
correlated state in the early stages with the initiation of 
a superradiant pulse from vacuum noise. Subsequently 
this is followed by a nonlinear regime where the super- 
radiant pulse is fully developed and its interaction with 
atoms becomes completely nonlinear. The decay of the 
pulse then occurs in a final linear dynamic stage. At 
the peak of SR, the collective radiation time of the sys- 
tem Tr = (STr/nA^L) T - IQ-^T - lO^i" s, becomes 
much smaller than the normal spontaneous emission time 
T ^ 60 ns for typical systems, where n is the density of 
atoms in the excited state, A is the resonant transition 
wavelength. 

For an elongated radiating sample, as the condensate 
along the z-axis being discussed here, superradiant emis- 
sion occurs dominantly along the ±z directions, i.e., 
emitted photons leaving the cigar shaped sample mainly 
from both ends. The corresponding spatial modes are 
called end-fire modes. They are perpendicular to the 
propagation direction of the pump laser beam. Due to 
momentum conservation for individual scattering events, 
the emission of end-fire photon is accompanied by col- 



lective recoils of condensate atoms. The momentum of 
recoiled atoms is significantly larger in magnitude than 
the residue momentum spread of the trapped condensate, 
thus gives rise to distinctly different islands of clouds 
upon on free expansion. These are the so-called con- 
densate side-modes. When the side-modes are occupied 
significantly, they serve as new sources for higher order 
SR, or sequential SR. They, too, emit end-fire mode pho- 
tons and contribute to the next order side-modes. The 
resulting pattern for atomic distribution after expansion 
as depicted in Fig. [U corresponds to what was observed 
for a certain choice of pump power and duration in the 
first BEC SR experiment 8] . The directions of the emit- 
ted end-fire mode photons and the corresponding re- 
coiled side-mode condensate bosons are indicated with 
the same line type. 



B. Effective Hamiltonian 

The effective second quantized Hamiltonian, govern- 
ing the dynamics of sequential SR system, is derived as 
follows. Due to the large energy scale difference between 
the center of mass (CM) dynamics for the atoms ('--^MHz) 
and the internal electronic degrees of freedom ('--^PHz), 
we can treat separately their respective motions. As in 
Ref. [2^ , the Hamiltonian of an atomic condensate with 
two level atoms interacting with a near-resonant laser 
pump, takes the following form 

H ^ j d?v^l{v) (-|^V2 + y,g(r)) ^,(r) 

+ j d?v4,l{v) (^-|^V2 + I4c(r) + hLo}j Mr) 
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%*(k)e-*-^t(r)4^,(r) + H.c.l ,(1) 



under the dipole and rotating wave approximations. We 
have further neglected the static atom-atom interac- 
tions. The first two terms are the atomic Hamilto- 
nians for the CM motion in their respective trapping 
potentials (T4g(r), Vto(r)) of the internal states. The 
atomic fields, described by annihilation (creation) oper- 
ator V^p,e(r) (ipl^ei^)), obey the usual bosonic algebra. 
Eg = and Eg = hA = h{LLjA — ^o) are the electronic 
excitation energies of the atom in the rotating frame de- 
fined by the pump laser field. The third term comes from 
the free electromagnetic field, while its interaction with 
the atoms is described with the last term (Haf), which 
includes both the laser photons and the scattered pho- 
tons. The operator Sk (0^) annihilates (creates) a pho- 
ton with wavevector k, polarization et, and frequency 
ujk = ck — ujQ (again in the rotating frame with frequency 

ujo). 5(k) = [c|k|dV2?ieo(27r)3]^/^ |k x x| is the dipole 

coupling coefRcient, with d = (e|r|(?) the matrix element 
for the atomic dipole transition. 

In typical SR experiments [8], the detuning A ~ 10^ Hz 
is much larger than both the CM motion energy scale 
~ 10^ Hz or the Rabi frequency ~ 10* Hz, the ex- 
cited state operator can be eliminated adiabatically via 
replacing V'e(r) ~ — (1/A) (/ dk5(k)e*''''"ak) V''s(r) in the 
equations of motion, yielding an effective Hamiltonian 



H = 



I d'r^lir) (^-|lv2 + ytg(r)) ^,(r) 

d^rdVk'~g{k, k', r)^J(r)a|^ak'V's(r), (2) 



with g(k, k', r) ~ g*{'k)g(k') exp (— «(k — k') • r), propor- 
tional to the effective coupling between the absorbed and 
subsequently emitted photons. 

The atomic field operators can be expanded in 
terms of the quasi-particle excitations of BEC V'g(r) = 
^q(q|r)cq, as described in Ref. with Cq (c^) anni- 
hilating (creating) a scattered boson in the momentum 
side-mode q in the form (r|q) = (/)o(r)e*''''. The ini- 
tial condensate mode is described by the spatial wave 
function 00 (r). The quasi modes for excitations approxi- 
mately form an orthonormal basis because (q|q') — Sq_qi. 
In the second quantized form within the side-mode rep- 
resentation Eq. ([2) becomes 

H ~ hujqcl^Cq + / d'^khuikol.ak 
q 

"K^/ d'kd^k'g*{k)g{k')pq.q,{k,k') X 



where pq,q/(k,k') = / dr|0o(r)Pe''^''+''^"^'''+'''^l '" is the 
structure form factor of the condensate density, which 
is responsible for the highly directional emission of the 
end-fire mode photons, ujq = ?i|qp/2m is the side-mode 
energy at the recoil momentum of q. The first two terms 
in Eq. ^ are diagonal in their respective Fock spaces 
and can be omitted by performing further rotating frame 



transformations Cq 



CqC *'^<i and dk 



Thus, the effective Hamiltonian takes the form 

^^-^E / ^^'kd3kV(k)g(k')Pq,q'(k,k')x 
q,q' 

X c^al^ak' Cq' e<"'^+"« (4) 

In a sufficiently elongated condensate, large off-axis 
Rayleigh scattering is suppressed with respect to the 
end-fire modes [31| . The angular distribution of the 
scattered light is sharply peaked at the axial directions 
(kg = ±kez), if the Fresnel number is larger than one 

= W^/LXq > 1 at the pump wavelength Aq for a con- 
densate of length L and width W [13] ■ This makes it 
possible to consider only the axial end-fire modes. To 
investigate sequential SR, we further take into account 
the first order side-modes at q = ko ± kg and the second 
order side-mode at q w 2ko. The rest of the side-modes 
are assumed to remain unpopulated [211] . The Hamilto- 
nian ^ that originally contains the contributions from 
all the side-modes and the end-fire modes as well as the 
laser field then reduces to the following simple model 



H 



doco + c^d^doco 



+C2dLdoc_ + C2d^doc+ j + H.c, (5) 



with g = g{ke). We have adopted a shorthand nota- 
tion where d± = d±k^, do = dko, c± = C(ko±kc)i and 
C2 = C2ko- This is the model Hamiltonian involving 
the interplay of the four atomic side-modes with three 
photonic modes. Before we further discuss and reveal 
the built-in entanglement swap mechanism for EPR-type 
quantum correlations in this model Hamiltonian, in the 
next section we shall briefiy review continuous variable 
entanglement and extend its criteria to our case. 



III. CRITERIA FOR CONTINUOUS VARIABLE 
ENTANGLEMENT 

The existence of continuous variable entanglement is 
determined by a sufficient condition on the inseparability 
of continuous variable states as given in Ref. [32] . If the 
density matrix of a quantum system is inseparable within 
two well defined modes [H, [s^, or the two modes are 
entangled, the total variance of EPR-like operators, u — 
\c\xi +X2/C and v = \c\pi — ■p2/c, satisfies the inequality 



Xc]jdj^dk'Cq', 



(3) 



(Ail' 



(Av') < 



(6) 
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for a real number c, where xi^2 = ('^1,2 + 2)/^^ 

Pi, 2 = ('^1,2 — o\ 2)/*V^ are analogous to position and 
momentum operators as in the case of a simple harmonic 
oscillator. The indices corresponds to mode numbers. 
Defining the inseparability parameter 



(7) 



the presence of continuous variable entanglement is then 
characterized by the sufficient condition A < 0. For the 
two modes to be entangled, it suffices to find only one 
value of c that leads to A < 0, and hence c can be taken 
at which A is minimum. 

The parameter A we adopt here clearly corresponds 
to an entanglement witness, but not an entanglement 
measure. Because the states for more negative A do not 
necessarily correspond to more entangled states. 

The total variance of the EPR operators are bounded 
below by the Heisenberg uncertainty relation (A-u^) + 
{Av"^) > ~ l/c^|- Thus, A has a lower bound A 



low 



l/c^ 



l/c2 



After minimization with respect to c, a more explicit 
expression of A can be given as 



A = 2 (c^(a{ai) + {(x\a2) /c^ + sign(c)(aia2 + a\a\) 



2/ 



where = 



(8) 



1/2 



_((aLa_)-|(a_)p)/((a;a+)-|(a+)n 

with sign(c) — —sign [Re {(a+a_)} — a+Q;_ + /3+/3_], 
a± — Re{(a±)} and /3± — lm{(a±)}. 

In the system we study here, the bosonic mode op- 
erators will be either end-fire mode pairs ai 2 = a± or 
end-fire modes and first side-modes, ai = a±,a2 — c^. 
Unlike other model investigations [33] of EPR-like corre- 
lations based upon A, we need to keep track of the (u)^ 
and (■i))^ terms, because (xi,2) and (^1,2) do not necessar- 
ily vanish for our model during time evolution. Further- 
more, since the time evolution of the two end-fire modes 
are symmetric in our case, we find — 1 and Aiow = —2. 

In the remainder of this paper, we examine the time 
evolution of the continuous variable entanglement wit- 
ness \{t) both for the opposite end-fire modes and for 
the end-fire modes with side-modes. This study is ex- 
pected to provide insight into the temporal development 
and the swap of quantum correlations between different 
sub-systems/modes. The following section is aimed at 
establishing an intuitive understanding of how EPR-like 
correlations between opposite end-fire modes are built 
up. 



IV. AN ENTANGLEMENT SWAP MECHANISM 



In sec. I VII we will exhibit the numerical results for the 
time evolution of the entanglement parameter A(t), gov- 
erned by the Hamiltonian ([5]). We will observe that, there 



exists regions in time where A becomes negative, i.e., con- 
clusive evidence for the presence of entanglement during 
dynamical evolution. In this section, we hope to provide 
an intuitive understanding to support the result revealed 
through the numerical approach. We will show that it 
is due to the presence of an inherent swap mechanism 
which leads to the generation of the EPR photon pair. 
We shall examine the dynamical behavior of the system 
in two different time regimes, the early times when the 
first side-modes just start to grow and the later times 
when the second order side-mode contributes to the dy- 
namics. 



A. Early Times 

In the initial stage, occupation of the second order 
side-mode (|c2)) can be neglected. During this initia- 
tion period of the short-time dynamics, the number of 
atoms in the zero-momentum state can be assumed un- 
depleted cq ~ -v/TVe**^^ with a constant N standing for 
the number of condensed atoms, like in the treatment 
of degenerate parametric processes. Since the pump is 
very strong, and the number of pump photons are much 
larger than the number of condensate atoms M ^ N, it 
can also be treated within the parametric pump approx- 
imation So « VMe*^° as undepleted. Thus, the initial 
behavior of the system is governed by the Hamiltonian 



Hi = -hxi 



H.c. 



(9) 



with xi = ^/NM\g\'^/A and 9i ^ 0^ + 0i is the ini- 
tial phases difference. This form of Hi is exactly the 
same as that of two uncoupled optical parametric ampli- 
fiers (OPAs). It allows for the growth of the first-order 
side-modes J30| as well as the entanglement of side-mode 
atoms with the end-fire mode photons [34| . The solution 
to Hi in the Heisenberg picture is given by the following 
time dependencies of operators |35l] 

a±{t) = cosh(xii)a± -I- ze'^sinh(xit)c^, (10) 
c±{t) = cosh(xii)c± -I- ie*^sinh(xit)a?p, (11) 

where the operators without time arguments are at the 
initial time. 

The side-modes and the end-fire modes are initially 
unoccupied |a+, a_, c+, c_) = |0, 0,0,0), or taken to be 
in their respective vacuum states. The time dependencies 
for the populations of the side-modes and end-fire modes 
come out as {I±) = (n±) = sinh^(xii), analogous to the 
classical results [35| . Evaluating correlations between the 
two end-fire modes (e), we find 



A = Acc = 4sinh^(xit), 



(12) 



which is always positive Aec > 0. On the other hand, 
the correlation between the end-fire mode (e) and side- 
mode (s), scattered in the opposite directions, takes the 
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following form 

Aso = 2 [2sinli=^(xit) - |sin(0i)|sinh(2xit)] , (13) 

which starts with Asc (t) = and evolves down to Ago (^o ^ 
2) ~ —2, the lowest possible value of Aiow = —2 imposed 
by the Heisenberg uncertainty, at | sin(6'i)| ~ 1. 

In Figs. [3Ka) andj^Jb), this is further supported by 
more elaborate results from numerical calculations of A 
and Ago. Same conclusions can be seen for the analytical 
results over the range t = 0-t = 0.3ms in Fig. ^a.) 
and (|13p over the range t = 0-0. 2ms in Fig. ^h). 



B. Later Times 

At later times, the first-order side-modes become ef- 
fectively populated, giving rise to noticeable second se- 
quence of SR from the edges of these side-mode conden- 
sates. In this case, the occupancy for the |co) mode is not 
an important issue, but the |c2) mode becomes populated 
due to the second order SR. 

We construct an approximate model by assuming that 
the occupation of |c2) is not changing too much, or ef- 
fective treating it as in in the steady state with £2 ~ 
y^/Y2g-«02^ iV2 is the number of atoms in the |c2) state. 
The later stage dynamics of the system, where the 2"^^^ or- 
der SR is effective, is then governed by the model Hamil- 
tonian 



-hx2 



H.c. 



(14) 



with X2 = y/N2M\g\^/A and 6*2 = 6*0 + 4>2- As before we 
again neglect the depletion of the pump oq ~ \/Me*^" . 

This model H2 is also exactly solvable. The time de- 
pendencies of the annihilation operators in the Heisen- 
berg picture are 



a±(t) = cos(x2At)a± 



^sin(x2At)c±, (15) 



c±{t) = cos(x2At)c± +ie-"^^sin(x2Ai)a±, (16) 

where t > t^, the operators without time arguments are 
a,t t = tf), and At = t — to- We can approximately con- 
nect these two models together into smooth temporal dy- 
namics if we use the solutions of Hi as the initial state 
for dynamics due to H2 so that a±(to) ^^nd c±(to) are 
calculated a,t t = to from the equations pO]) and (fTTj) . 
respectively. We define as the time at which all the 
I Co) atoms are scattered into the side-modes, and thus it 
is determined from sinh^(xiio) = A^/2. 

In this later dynamical stage, the entanglement witness 
parameter in between the end-fire modes (e) is evaluated 
to be (A = Acc) 

X{t) = 4sinh2(xiio) 

- |cos(0)sin(2x2At)|sinh(2xiio), (17) 

where 9 = 9i + 62- When | cos(0)| ~ 1, A evolves from 
2N down to the minimum possible negative value of 



Alow = —2, at At = TV / 4x2- An analogous calculation for 
entanglement between the end-fire mode (e) and side- 
mode (s) gives 



Xsc{t) = 4sinh2(xito 



2|sin(^?2)cos(2x2At)|sinh(2xiio), (18) 

which starts at Asc(to) — —2 and increases to values of or- 
der ~ N, for proper choices of 92- Many of these features 
revealed from simple analytic models find their parallels 
in numerical solutions as displayed in Fig. [31 

The results from the two model Hamiltonians are found 
to depend on the initial phase difference between 9i and 
6*2, but not the individual phases. Such a phase depen- 
dence of the results is analogous to the cases of paramet- 
ric down conversion and the two-mode squeezing [35| . 
The phases introduced in the second stage refiects the 
accumulating temporal phase difference of the operators 
through the time evolution. In the numerical calculation 
it is sufficient to assign initial phases for the pump laser 
and the condensate or just their difference. 

Without any detailed analysis, simply consider the be- 
haviors of Eqs. and ([TT]) instead, one can already 
appreciate the built-in entanglement swap mechanism 
within the superradiant BEC in action. The entangle- 
ment created between the side-mode and end-fire mode 
Eq. (|13p in the initial stage, is swapped to entangle- 
ment between the two end-fire modes Eq. (fT7)) . due to 
the 2"'^ order SR. The model Hamiltonian Hi couples 
the \a±) <-> |c^) modes, but leaves |a+) ^ |a_) modes 
decoupled at the initial times. The model Hamiltonian 
H2, at later times, couples the \a±) ^ \c±) states. Two 
noninteracting modes \a+) <-> |a_) are coupled through 
their common interaction with the same side-mode, and 
become entangled due to the swap mechanism. 



V. NUMERICAL CALCULATION OF THE 
ENTANGLEMENT PARAMETER 

We study the dynamics of the entanglement parame- 
ter A(t) and the accompanying populations for the fields 
(/o(t),/±(t)) and the atomic states (no,n±(t),n2). Their 
complete temporal evolution is governed by the Hamil- 
tonian Eq. Our calculation will be numerically ob- 
tained, aided by a decorrelation approximation that ne- 
glects higher order correlations. The numerical results 
will be illustrated and discussed in the next section. 

The entanglement parameter A, given in Eq. ([5]), is 
determined by the expectation values of both a± oper- 
ators and their products. Their equations of motion in 
operator forms can be derived from the full Hamiltonian 
Eq. ^ . The dynamics of two operator products is found 
to depend on four operator products; and four operator 
products depend on six operator products, and so on 
so forth. Such a hierarchy of operator equations is im- 
possible to manage in general. We therefore resort to a 
decorrelation approximation that truncates it to a closed 



6 



form. The usual treatment of this kind [36!| for the SR 
system closes the chain early by a simple decorrelation of 
atomic and optical operators, which is clearly inappro- 
priate when entanglement swap is to be studied. 

We adopt a decorrelation rule that factorizes conden- 
sate and the second order side-mode operators in oper- 
ator products. Since quantum correlations between the 
condensate and other modes are expected to be weak 
due to the almost classical, coherent-state-like nature for 
the condensate and its diminishing population when the 
second order side-mode is significantly populated at later 
stages of dynamics. Operators for the pump photons will 
also be factorized, again relying on the almost classical, 
coherent state nature of the pump field. 

Our approach makes it possible to keep quantum cor- 
relations between the end-fire modes and the interme- 
diate side-modes. The hierarchy of equations is closed 
under {xyz) « {xy){z), with x,y G {1, c±, clj_, a±, alj_} 
and z G {co, Cq, C2, C2}. The resulting equations, govern- 
ing the dynamics of the expectations, are given in the 
Appendix. These equations are solved numerically. 

For the initial conditions, both the end- fire modes and 
side-modes are taken to be their vacuum Fock states 
while the laser and the condensate are in coherent states. 
We consider a system with typical parameters of a con- 
densate with number of atoms = 8 x 10^ and a pump 
with M = 2 X 10* photons. Additionally, phenomenolog- 
ical decoherence rates arc introduced by assuming the 
same damping rates [3ll | for the atomic and photonic 
modes. The decay rates are obtained from the effective 
decay of the experimentally measured contrast for the 
atomic density distribution pattern 3. In addition, we 
also explored an interesting scheme where the coupled 
equations (|AHIA18|) were solved, in the presence of phe- 
nomenological damping, for an initial two-mode squeezed 
vacuum (for the end-fire modes) with a squeezing param- 
eter ^ = r exp (iOr). 

VI. RESULTS AND DISCUSSION 

In sec. llVl we discussed the origins of the entanglement 
swap in sequential SR. In this section, in order to provide 
for a more detailed and quantitative understanding, we 
present results obtained from numerical calculations. We 
will discuss the time evolution of the entanglement pa- 
rameter X(t), between the two end-fire modes, within the 
parameter regime of the experiment Q. At first, we will 
disregard decoherence and examine the nature of fully co- 
herent sequential dynamics. We will show that A attains 
negative values, confirming the presence of entanglement 
due to the swap mechanism as we intuitively discussed in 
the previous section. We then introduce effective damp- 
ing rates specific to the experimental situation. Finally, 
we will examine the effect of initializing the quantum dy- 
namics of our model system from in a two-mode (end-fire 
modes) squeezed vacuum, in the presence of decoherence 
and dissipations. We will end with investigations of the 



dependence of correlations on the number of condensate 
atoms. 



A. Dynamics of Entanglement 

In Fig. [51 we plot the temporal evolution of optical 
field intensities and atomic side-mode populations. The 
plot is found to be totally symmetric with respect to 
t — tc — 0.35 ms. The peak in the intensity after tc is 
the analog of the Chiao ringing [s^l • In the experiments 
such a complete ringing cannot be observed due to the 
presence of decoherence, which will be taken into account 
in our subsequent calculations. 

In Fig. [31 we plot the temporal evolution of entan- 
glement parameters Aoc (photon-photon) and Xgc (atom- 
photon) over the population dynamics, depicted in Fig. 
[21 The lower panel of Fig. [31 demonstrates the swap dy- 
namics. The initial atom-photon entanglement (Ago) is 
seen to evolve continuously into entanglement between 
the two end-fire modes (A). Both the parameters Age and 
A is found to be able to reach down to the lowest possi- 
ble value, Alow = —2, set by the Heisenberg uncertainty 
principle (as in sec. IIIip . The complete numerical results 
match well with the analytic solutions, discussed previ- 
ously in Sec. IIVI for the model Hamiltonians Q at early 
times and (fH]) at later times. 
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FIG. 2: (Color online) The temporal evolutions for atomic 
side-mode populations and optical field intensities. I±, n±, 
no, and 722 denote occupancy numbers of bosonic modes \a±), 
c±), I Co), and |c2), respectively. n±{t) and I±{t) overlap 
except for a short time interval near t — tc = 0.35ms. Notice 
that no and n2 are scaled for visual clarity. 

In the time interval of t = 0-0.30 ms, we find the sys- 
tem is dominated by the 1*^' sequence of SR. The atomic 
condensate, initially in the zero- momentum state |co), 
is pumped into the 1*** order side-modes \c±). This is 
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FIG. 3: (Color online) (a) The temporal evolutions of atom- 
photon (|a±) ^ |ct)) ^i^^ photon-photon (|a+) « 
mode correlations as evidenced by the entanglement parame- 
ters Asc and A = Acc, respectively. Accompanying population 
dynamics (Fig. [Sjl is also plotted in the background, (b) an 
expanded view of the early time dynamics for Ase and A; (c) 
an expanded view of A around tc — 0.35 ms. 



the reason for the overlap of n±{t) with I±{t) during 
this interval. Due to the interaction between the side- 
modes and the end-fire modes, scattered into opposite 
directions, Age becomes negative in this region. 

When the \c±) side-modes become maximally occu- 
pied at about t = 0.30 ms, the 1^' sequence of SR is com- 
pleted. At this time, these side-modes are sufficiently 
populated to give rise to the 2"^^ sequence of SR. In the 
interval t = 0.30-0.35 ms atoms in the side-modes \c±) 
are pumped into the 2"'^ order side-mode |c2). The ma- 
jority of the populations, however, oscillates back to the 
I Co) mode because of the more dominant Rabi oscillation 
between the |co) and \c±) modes. Two other reasons also 
contribute to the repopulation of the condensate mode: 
first, the neglect of the propagational induced departure 
of the end-fire mode photons from the atomic medium; 
and second, the neglect of the other two 2"'^ order side- 
modes |c2ko±2ke) for atoms to get into. Two end-fire 
modes get indirectly coupled by the entanglement swap- 
ping and between t — 0.30-0.35 ms A(t) gradually be- 
comes negative. 

Entanglement of the end-fire modes arises a,t t — tc — 
0.35 ms, when the |c2) mode is maximally occupied as 
shown in Fig. [3) The minimum value of A, which occurs 



at t = tc, is found to coincide with the maximum value 
of n2{t). 

When t > tc, however, due to our limited mode ap- 
proximation of not including even higher side-modes, we 
cannot study any effects which could potential give rise to 
higher order correlations, such as the onset of the 3'''^ se- 
quence of SR. The oscillatory Chiao type ringing revivals 
in the present result after t > tc mainly arise from the ex- 
clusion of decoherence, dephasing, dissipations, and the 
higher order side-modes in the model system. In the 
present work, we limited ourselves to a particular side- 
mode pattern as actually observed in available experi- 
ments 0. Despite its simplicity, we find our model can 
reasonably explain effects of decoherence and dephasing 
on the entanglement dynamics, which is further illus- 
trated in the next subsection. 



B. Vacuum squeezing and Decoherence 

The introduction of experimentally reported decoher- 
ence rate of 7/27r = 1.3 x lO** Hz phenologically into the 
dynamical equations for the coupled system is found to 
not change the nature of the entanglement and swap dy- 
namics significantly, which is supported by the numerical 
results shown in Fig. [D We find that A can still become 
negative in certain time window, although it now stops 
short of reaching the theoretical lower bound of —2. 

In the lower panel in Fig. 3] the temporal window for 
the negative values of A, or the presence of entanglement 
is found to become narrower and the minimum value of 
A, Amin, is now somewhat larger for stronger decoher- 
ence, as may be expected. According to sec. IIIIl a less 
negative value of A does not necessarily imply less entan- 
glement, because A is simply an entanglement witness pa- 
rameter but not an entanglement measure. On the other 
hand, it is still beneficial to aim for lower values of A, 
because the numerical results we obtain associate lower 
values with longer entanglement durations, and further- 
more more tolerant to decoherence, which means photon- 
photon entanglement can withstand the higher decoher- 
ence rates. 

For this aim, we choose to consider end-fire modes 
which are initially in two- mode squeezed vacuum states. 
The lower panel in Fig. 2] shows that an initially two- 
mode squeezed vacuum, for the end-fire modes, can in- 
deed compensate to a certain degree for decoherence. 
This shows that, initially induced two-mode squeezing 
(or entanglement) in between the end-fire modes en- 
hances their subsequent entanglement after the entan- 
glement swap. 

This observation can be interpreted as follows based 
on the numerical results. Any initial correlation between 
the end-fire modes is lost in the early dynamical stage 
where the end-fire modes are entangled with the first 
side-modes. The presence of initial correlation, how- 
ever, causes the resultant atom-photon entanglement to 
be more resistant to decoherence. As a result, photon- 
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FIG. 4: (Color online) (a) The temporal evolution of atomic 
and field mode populations, and of entanglement parameters. 
A decoherence rate of 7o/27r = 1.3 x 10*H[z is introduced 
without any initial squeezing, (b) An expanded view of the 
dependence of A on decoherence rate 7 and squeezing param- 
eter r around tc = 0.38 ms. 



photon correlations established by swapping from the 
atom-photon correlations in the subsequent dynamical 
stage also become more resistant to decoherence. 

Finally, we examine the influence on A from the num- 
ber of atoms in a condensate. We find that, as illustrated 
in Fig. [5l Amin becomes more negative for larger conden- 
sates. In the small condensate limit, Amin is found to 
decrease linearly with N when the Fock vacuum is con- 
sidered as initial conditions for other modes. The lower 
limit of —2 is never attained. When a small amount of 
initial squeezing is introduced, however, A can be brought 
down to theoretical minimum of —2. It approaches —2 in 
the large condensate limit with or without any help from 
initial squeezing in the two end-fire modes. 

In addition to the amplitude of squeezing parameter, 
its phase could also influence Amin- In Fig- El we plot 
the minimum value of the entanglement parameter as 
a function of the phase and amplitude of the squeezing 
parameter ^ = re**''. We performed this study for a 
small condensate with = 100 atoms and ignored the 
phenological decoherence. We find that the most efficient 
enhancement occurs along the line = t^- For larger 
condensates we find that the center of Fig. O where 
A > 0, spreads out to the 0^ = and Or = edges, as 
N is increased. Entanglement is enhanced mainly along 
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FIG. 5: The dependence of Amin on A'' in different scales. Solid 
lines are for an initial coherent vacuum (r = 0) and dashed 
lines are for a squeezed vacuum (r = 0.005 and 9 = n). 
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FIG. 6: (Color online) The dependence of Amin on r and Or 
for N = 100. Amin shows a mirror symmetry for 9r > tt. 



VII. CONCLUSIONS 

We investigate photon-photon entanglement between 
the counter-propagating end-fire modes, of a sequentially 
superradiant atomic Bose-Einstein condensate. We cal- 
culate the temporal evolution of the continuous variable 
entanglement witness parameter for suitable realistic ex- 
perimental parameters in the cw-pump laser regime 0], 
and find that EPR like correlations can be generated be- 
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tween the oppositely directed end-fire modes, despite the 
fact that they do not directly interact. 

The generation of entanglement is shown to be due to 
a built-in entanglement swap mechanism we uncover in 
the sequential SR system. It is shown that end-fire mode 
photons become entangled immediately after the second 
sequence of the superradiance. In the second sequence, 
one of the end-fire modes interacts with the side-mode, 
with which the other end-fire mode has already interacted 
before in the first sequence. This mechanism allows for 
swapping the entanglement established between the end- 
fire modes and the side-modes in the first sequence, to 
the entanglement of the end-fire modes per se. 

Increasing the number of atoms in the condensate, or 
initializing superradiance with a two-mode squeezed vac- 
uum (for the end-fire modes) , are found to be beneficial to 
the efficient construction of entanglement between end- 
fire modes, via the increasing of entanglement durations 
and making the entanglement more tolerant to decoher- 
ence. 

The initial phase difference of the incoming pump laser 
and the condensate, the phase and the amplitude of the 
squeezing parameter for the end-fire mode vacuum, the 
number of atoms in a condensate and its geometric shape, 
all play certain roles in order to achieve the optimum 
ERP-like correlations in between the end-fire modes and 
these parameters are discussed in detail in the present 
manuscript for the cases of both small and large conden- 
sates. 
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APPENDIX A 

We calculate temporal evolution of entanglement pa- 
rameter X{t), given in ([5]), starting from the Heisenberg 
operator equations, obtained from (O. We evaluate the 
expectations for both single operators and two opera- 
tor products. We arrive at a closed set from the expec- 
tations via performing decorrelation approximation, in 
parallel with the development and understanding of the 
swap mechanism (Sec. IIVI) . 

The resulting closed set of equations for expectation 
values are given through Eqs. (|AmA18| . where time is 
scaled by frequency a = g'^ /2A with 5 ~ 2 x 10'^ Hz, while 
operators are not scaled, a is related to 7 of Ref. [31[, as 
7 = \[Ma = 10.7 Hz. Phenological decay rates can be 
introduced in Eqs. (|AHIA18|) by scahng 7^ = 1.3 x lO'^Hz 
with Of. However, since the decay rates are introduced, in 
[sH . for three-operator products, we use 7_l/3 for single 
operators and 27^/8 for two operator products. We have 
also checked the parallelism of our density dynamics with 



3l( . which are in good agreement with the experiment 
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dt 
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dt 



(Al) 



(A2) 
(A3) 
(A4) 
(A5) 

(A6) 



i(ao)((4)(co) + (cJ)(c±)j, 
z((a_c+)(4) + (a_cL)(c2) 
(a+C-)(4) + (a+c+>(c2)), 
*((«o)(a^p)(co) + (ao)*(a±)(c2) 
*(ao)*((a-c+) + (a+c_)^, 
i{ao)[{alc^) + (a^c+>), 

i(ao) (^(a_c^)(co) + (a_c+)(c^) 
(a+c;>(co> + (a+c_>(4)), 

i((ao)(aLc^)(co) + {ap) {a^_c+) {c\) 
(ao)*(a+c+)(4) + (ao)*(a+cL)(c2)),(A7) 

2*(ao) ((a±c^) (cq) + (a±c±) (c^ )) , (A8) 

*((«o)(a-C_)(co) + (ao)*(a+c_)(c2) 
(ao>(a^c+)(co) + (ao)*(a_c+> (ca)) , (A9) 

i((ao)(aLc^_)(co) + (ao)*(a+c^_) (ca) 
(ao)*(a+c+>(4> - (ao)(a^c+)(4)),(A10) 

2i(^(ao>(a5pC±)(co) 

(ao>*(a±c±)(c2)), (All) 
*((ao)(4c^)(co) + (ao)(c±Czp)(cT) 
(ao)(a±a±>(co) + (ao)* (a±a^) (02)) (A12) 

«((«o)(c±c?p)(co) + (ao)(c±clt)(cJ) 
(ao>*(a±aT>(4> - (ao>*(a±4>(4>)A13) 

i((ao)(c±4)(co) + (ao)(4)(4> 
(ao>(a±4>(co> + (ao)*(4>(c2)), (A14) 



iy - (ao)*(aTC±>(4> - (ao) (a±c±) (4) 
(ao)(a^4)(co) + (ao)*(a±cL)(c2)),(Al5) 



(4co) 
dt 



i[ - (ao)(aLc^)(co) - (ao) (a^cL) (cq) 



t J 
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+ (oo)*(a_c+)(4) + (oo)*(a+c_)(4)VA16) (a±a±) ./ , \/ t\ / \*/ t\/ \ 

= - (ao)*(a-cL)(c2) - (ao)*(a+c^)(c2) + (ao)(aLc^)(co> + (ao)(aLc±)(4)) . (A18) 

+ (ao)(alc_)(4) + (ao)(a;c+)(4)), (A17) 
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